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ATOMIC OXYGEN EFFECTS ON METALS 


by 


Albert T. Fromhold 
Professor of Physics 
Auburn University, AL 36849 


ABSTRACT 

The present work addresses the effect of specimen 
geometry on the attack of metals by atomic oxygen. This 
is done by extending the coupled-currents approach in 
metal oxidation to spherical and cylindrical geometries. 
Kinetic laws are derived for the rates of oxidation of 
samples having these geometries. It is found that the 
burn-up time for spherical particles of a given diameter 
can be as much as a factor of 3 shorter than the time 
required to completely oxidize a planar sample of the 
same thickness. 
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INTRODUCTION 


The failure of metal structures often occurs at 
cracks, corners, and crevices preferentially to planar 
surfaces. The present work addresses these geometric 
effects of the attack of atomic oxygen on metal 
surfaces. 


OBJECTIVES 

The objective of this work was to extend the 
coupled-currents approach for the oxidation of planar 
material surfaces to non-planar surfaces, specifically, 
surfaces having spherical and cylindrical symmetry. 


CONCLUSIONS AND RECOMMENDATIONS 

1. The burn-up rate can be significantly increased for 
non-planar geometries relative to planar geometries. 

2 . The theoretical development should be extended 
to include stress effects. 
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I. Concept of Steady-State Oxide Growth for Non-Planar Geometries 

The focal point of this presentation is the effect of material 
sample geometry on the rate of oxidation. Because the oxidation rate 
is dependent upon the atom or ion current to the location where metal 
oxidation is occurring, it is necessary to evaluate the effects of 
sample geometry upon such currents. Frequently oxide formation occurs 
under fixed experimental conditions of ambient temperature and oxygen 
pressure, and in such situations the currents in the system are usually 
nearly steady- state . 

The steady- state is a condition in which no change in the local 
value of the current takes place with time. Though thermodynamic 
equilibrium is a special case of the steady-state, the current then 
being uniformly zero, in general the steady-state differs from 
equilibrium. 

To elaborate, current flows in response to driving forces, so 
currents are to be expected so long as the net driving force differs 
from zero. Currents may lead to a change in one or more components of 
the local driving force, in which case the net driving force may 
approach zero in time. Then the net current will approach zero in 
time, and a state of equilibrium is thereby attained. This evolution 
of the system can be expected for many, but not all, situations. A 
specific counter example is that of an externally-applied driving 
force and a closed- loop circuit around which the current can be 
sustained practically indefinitely, so the equilibrium state is not 
approached with time. In fact, there will be initial transients when 
the current is established, and it is the steady-state which is 
approached with time. Thus for some systems the long-time asymptotic 
limit on a laboratory time scale can be the steady-state, in contrast 
to equilibrium. 

In the steady- state limit, the actual values of the local currents 
will depend directly upon the local values of the conductivity, as well 
as upon the local driving forces. The critical distinction between the 
steady- state and equilibrium, however, is in actuality one of net local 
driving force rather than actual value of the local current, even 
though in practice it may well be the measurable current that is used 
as the threshold criterion for equilibrium. Theoretically the values 
of the net local driving force must be zero everywhere in the system at 
equilibrium, and this naturally leads to the condition of zero current 
everywhere in the system, regardless of the value of the local 
conductivity . 

The space -dependence of the current is likewise pertinent. For 
non- steady- state situations, the current can certainly vary locally 
with position. In the steady- state, however, the .situation is not 
quite so intuitively clear. The local driving force or the local 
conductivity may vary with position initially, but the local driving 
force may adjust to compensate for the spatial variation of the 
conductivity. Care must be taken in the analysis of such cases. The 
critical aspect of the problem is that the spatial variation of the 
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current in the limit of the steady- state must be such that neither 
electric charge nor diffusing species concentrations build up or 
deplete anywhere within the system as the flow is maintained. This 
requires that the net current flowing into any region of space be equal 
to the net current flowing out of that region of space. Thus there 
must be no time -dependence of the local concentration of the mobile 
species. In vector terminology, the divergence of the current J in 
question must be zero , 

(1) -a c/at = V-J = 0 , 

where C is the local mobile species concentration at time t. This 
condition will be shown to lead to the conclusion that even in the 
steady-state, the current can vary locally with position, depending on 
the geometry in question. Our attention is directed entirely to cases 
of thermal oxidation for which Eq. (1) represents a reasonably good 
approximation . 


II. Ionic Current and Oxide Growth Rate 


If the steady- state ionic current density can be evaluated for 
the geometry in question, then the value of that current at the 
location where the chemical reaction of metal oxidation is occurring 
(usually one of the oxide interfaces) will enable a determination of 
the oxidation rate. Specifically, the ionic particle current density 
J. at the reaction interface leads to an increase in the oxide film 
thickness L at that interface. The local oxide thickness formation 
rate at that location can be written as 


< 2 > ft - R i J i 


where R. is the volume of oxide formed per transported particle of 
ionic species labeled i. The volume V oxi( j e of oxide being formed 
because of this oxide thickness increase is given by the product of the 
thickness increase with the area, summed over the reaction zone where 
oxide film formation is occurring. Assuming this to be either the 
interface between the sample and the oxide (designated the metal- oxide 
interface) or else the interface between the oxide and the oxygen being 
utilized for the reaction (designated the oxide-oxygen interface ) , then 
the oxide volume formation rate is given by 


(3) 


d V o ilde 
dt 



J • dA . 

Reaction 

‘'Interface 


Note that the surface integral involves the vector dot product between 
the current density vector and the differential area vector dA on 
the reaction interface. The current density vector is parallel to the 
motion of the species, and the differential area vector is locally 
perpendicular to the surface area where the current density vector is 
being evaluated. The effect of taking the vector dot product is to 
obtain the product of the current density flow perpendicular to the 
surface area with that surface area in order to obtain the total 
particle flow through the surface. 
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III. Sample Geometry 


A. Domain of Application 

Since the oxidation rate of materials is dependent upon the ionic 
currents through the protective oxide covering the material, it is 
necessary to evaluate the effects of sample geometry upon such 
currents. Our attention is now directed to the spatial variation of 
the steady- state current for cases of spherical and cylindrical 
symmetry. Spherical symmetry should well characterize the oxidation of 
spherical particles and should yield an approach to the oxidation 
within spherical cavities in materials. Cylindrical symmetry should 
well represent the oxidation of wires, and should model reasonably well 
the oxidation of edges formed by intersection of two planar surfaces. 


B . Spherically -Symmetric Case 

Consider a uniform radial flow of current density J s of species s 
between two concentric metal spheres having radii a and b, with a < b, 
which can be called the inner and outer electrodes. If none of the 
species concentrations C s change at any vector position r with time t 
the condition necessary for the steady-state, then the total current i c 
of species s given by s 


(M I s = 


J (r) 

4ir e e rt IV ot s 

Radius r 


dA *= J g (r) [ 4?rr 2 ] 


must be the same when evaluated over any imagined sphere concentric 
with the spherical electrodes, assuming a < r < b. The spherical 
symmetry insures that J s is independent of the position of the area 
element dA on the sphere of radius r. That is, there is no preferred 
radial direction of flow in a spherically- symmetric situation such as 
this, so the current density is uniform over the surfaces of symmetry 
consistent with the physical system. This allows the ready evaluation 
of the surface integral in Eq. (4). 


The fact that the total current I s is independent of the radius r 
of the chosen sphere in Eq. (4) means that for a spherical shell volume 
element bounded by two such spherical surfaces, the current out of the 
volume element will equal the current into the volume element. Thus 
there will be no build-up or depletion of the concentration of the 
species within the volume element, as required by the steady- state 
condition. 


The conclusion to be derived from Eq. (4) is that the steady-state 
radial current density for the case of spherical symmetry varies 
inversely with the square of the distance r, 

(5) J g (r) « 1/r 2 . 

C. Cylindrically- Symmetric Case 

If the two metal electrodes are concentric cylinders instead of 
being concentric spheres, then the current will likewise be radial if 
the medium between the two electrodes is uniform. Equation (4) is then 
replaced by 
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(6) I s - 


ISymm e try 
Cylinder 
Radius r 


J g (r) dA - J g (r) [2jrrh] 

O f 


where h is the length of the concentric cylindrical electrodes having 
radii a and b (a < r < b) . The steady-state in this case likewise 
requires that I s be independent of r, where in this case r is the 
radius of any imagined cylinder of symmetry with the electrodes 
(a < r < b) . The symbol dA in this case is the area element on the 
cylinder of radius r. The uniformity of the current density J s , which 
because of the electrode symmetry is in the radial direction and has 
the same value at any point on a given symmetry cylinder, leads to the 
ready evaluation of the surface integral as shown in Eq. (6). There is 
no need to consider the flat ends of the imagined cylinder in Eq. (6), 
since there would be no contribution to I s from a surface integral over 
those regions because there is no current perpendicular to the flat 
ends of the cylinder. The separation distance h between the ends of 
the imagined cylinder does enter into Eq. (6) because the surface area 
of the curved portion of the cylindrical surface through which the 
current flows depends upon the length of the cylinder. 


The independence of the total current I s on radius r of the chosen 
cylinder in Eq. (6) means that for a cylindrical shell volume element 
bounded by two such cylindrical surfaces, the current out of the volume 
element will be equal to the current into the volume element. Thus 
there will be no build-up or depletion of the concentration of the 
species within the volume element, as required by the steady- state 
condition. 


The conclusion to be derived from Eq. (6) is that the steady-state 
radial current density for the case of cylindrical symmetry varies 
inversely with the distance r, 


(7) J g (r) cc 1/r . 


D . Planar- Symmetric Case 

For completeness, let us carry out the above analysis for two 
metal electrodes which are symmetrically- located parallel planes 
separated by some distance L. The current will be perpendicular to the 
electrodes if the medium between the two electrodes is uniform. 

Equation (4) is then replaced by 


(8) I s 


J (x) dA = J (x) A , 

S y mm etry S S 

‘'Rectangle at 
Position x 


where A is the total area of the rectangle. The steady- state in this 
case also requires that I s be independent of x, where x is the distance 
to the rectangle, as measured perpendicularly from one of the 
electrodes . 


The symbol dA in this case is the area element on the rectangle. 
The uniformity of the current density J s , which because of the 
electrode symmetry is in the perpendicular direction and has the same 
value at any point on a given symmetry rectangle, leads to the easy 
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evaluation of the surface integral as shown in Eq. (8). There is no 
dependence of I s on the position x of the rectangle, as noted from Eq . 
(6). In addition, there is no dependence on the electrode separation 
distance L. 


The independence of the total current I s on position x of the 
rectangle in Eq. (8) means that for a parallelepiped volume element 
bounded by two such rectangular surfaces, the current out of the volume 
element will be equal to the current into the volume element. There 
would be no current through the four ends of the parallelepiped because 
the current density is parallel to those planar ends. Thus there will 
be no build-up or depletion of the concentration of the species within 
the volume element, as required by the steady-state condition. 

The conclusion to be derived from Eq. (8) is that the steady-state 
current density for the case of planar symmetry does not vary with the 
position x in the medium, 

(9) Jg <x Constant , independent of poxition x . 

To summarize, the steady-state current density is independent of 
position for planar geometry, it decreases inversely as the radial 
position r for cylindrical symmetry, and it decreases inversely as the 
square of the radial distance for spherical symmetry. 

E. More Formal Treatment 

We can treat the problems of spherical, cylindrical, and planar 
symmetry more formally by introducing the corresponding vector forms 
for the divergence operator in the appropriate orthogonal coordinate 
systems . 


The divergence operator in the spherical polar coordinate system 
is given by 


(10) V-J = 

r z sin 8 

where the vector 


[sin » fj (r 2 J r ) + r fj (sin <J,> ] 


(11) J-rJ + 0 J. + $ J , 

r v <p 

in spherical polar goojdinateg r, 0, and <t> has the components J r , 

Jfl , and , with r, 0, and <f> being the unit vectors in that 
coordinate system. Because in our application, J represents the 
current which is entirely radial and moreover is a function of r only, 
J r = J r (r), Jg = 0, and - 0. Equation (10) for the divergence 
reduces in this case to 


<12 > V ' J - h <' 2 J r> • 

Applying the steady-state condition in the form of Eq. (1) then leads to 

<13) fcp V r >]- 0 ■ 

which integrates to yield 
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(14) r 2 J r (r)=X sph 

with X sph being the constant of integration. Thus we obtain J oc l/r^ 
for this case, in accordance with Eq. (5) previously deduced. By 
comparing this result with Eq. (4) we see that X sph =I s /47r . 


The divergence operator in the 
given by 

(15) V-J - i d -(r J r > + i 
where the vector 

(16) J = r J + 0 3 . + z J 

r u z 

in cylindrical polar c§ordina£es r, 6 , and z has the components J r , 

Jfl , and J z , with r, 0, and z being the unit vectors in that 
coordinate system. However, in cylindrical symmetry cases, the current 
density vector is radial. Therefore, in cylindrical polar coordinates, 
the component J r is the only non- zero component of the vector; the 
components Jq and J z are zero. Moreover, J r is a function of r only. 
Thus we can write 

(17) V-J - ifj(rJ r ). 

Applying the steady- state condition in the form of Eq. (1) then leads to 

< 18 > V r) ] - 0 ■ 

which integrates to give 


cylindrical polar coordinates is 


+ ^ (J Z> • 


(19) r J r (r) - ^ yl , 

with J£ yl being the constant of integration. Thus we obtain J oc 1/r 
for this case, in accordance with Eq. (7) previously deduced. By 
comparing this result with Eq. (6) we see that X cy 1 -I s /27rh. 


For completeness, let us also carry out the above analysis for A 
planar symmetry, where any variation is considered to occur in the x- 
direction, with the system being entirely uniform in the y and z- 
directions. The divergence operator in rectangular cartesian 
coordinates is given by 


( 20 ) 


V-J = 


dJ ? 

dx 



+ 


aJz 

dz 


where the vector 


(21) J = X J x + y j y + Z J z 

in rectangular cartesian^coordinates x, y, and z has the components 
J x , J y , and J z , with x, y, and z being the unit vectors in that 
coordinate system. However, in the presently-considered case, the 
current density vector is in the x direction and moreover is at most a 
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function of the x- coordinate . Thus 

aj x ( X ) 


(22) V-J 


3x 


Applying the steady-state condition in the form of Eq. (1) then leads to 
3J„(x) 

(23) ~ir~ - 0 • 

which integrates to yield 


(24) J x (x) - X plane , 

with K plane being the constant of integration. Thus we obtain J x 
independent of x for this case, in accordance with Eq. (9) previously 
deduced. 


In the following section we show how to proceed in obtaining the 
values of the total current. This is tantamount to evaluating the 
integration constants obtained above. 


IV. Driving Forces and Currents 


The general vector relation for the current density J s for species 
s in terms of the electrochemical potential u s for species s is 
determined by the vector gradient, 


(25) J a - C Vu , 
s s s 


where C g is the concentration of species s. The form for the gradient 
operator V is specific to the symmetry of the physical problem to be 
treated. It can be expressed, for example, in rectangular, spherical, 
and cylindrical coordinates, or in a more general form which applies to 
any orthogonal curvilinear coordinate system. In spherical polar 
coordinates, 


(26) Vu 

s 


a 3u s A i 3u s a 1 3u_ 

r — + d ~r~~ + 

dr r 38 v r sin 6 3<f> ’ 


while in cylindrical polar coordinates, 


(27) 


Vu = r 
s 


duo A 1 

t - 5 + e - 

dr r 


3u £ 

w 


+ z 


3 u £ 
3 z 


Denoting the proportionality factor for Eq. (25) by 0 , and 

using the standard form for the electrochemical potential, 5 

(28) u = u° + k T In C + q V , 
s s B s M s 

where u° is the reference state value of the chemical potential u s 
kg is Boltzmann's constant, T is the absolute temperature, q s is the 
charge per particle of species s, and V is the electrostatic potential, 
Eq. (25) becomes 


(29) J s 


0 C Vu = 
s s s 


k„T 0 VC 
B s s 


0 E C , 
s s ’ 
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where E is the electric field given by 

(30) E - - W . 

Comparison with the usual linear diffusion equation 

(31) J s - • D s VC s + ^ E C s , 

shows term-by-term agreement, with the diffusion coefficient being given by 

(32) D s = k B T » 

and the mobility being given by 

(33) „ s - q s « • 

The ratio of the two coefficients gives the Einstein relation 

(34) „ s /D s - q s /k B T. 

The electric field is of course zero in the above equations if the 
diffusing species giving rise to the oxidation process are uncharged. 

With this choice, the equations in our development for the rate of 
oxidation yield also the growth rate for non-planar geometries for the 
diffusion of uncharged particles. In this limj.t the results can be 
compared to those deduced by Wilson and Marcus . Generally the 
diffusing species are charged, however, in which case the full 
treatment is required. 


V. Electric Fields 

A. Generals Relations and Zero -Space -Charge Limit 

For charged particle diffusion, the driving force of the electric 
field is as important as the concentration gradient, both being ^ 
included as part of the gradient in the electrochemical potential. For 
cases of high symmetry, Gauss' law provides the easiest approach to 
finding the value and position- dependence of the field. A closed 
imaginary surface having the symmetry of the physical problem is used 
for carrying out a surface integral of the electric field, 


(35) 


c E * dA 

Imaginary 

‘'Surface 


Q 


where Q is the net charge within the closed surface that is chosen 
for the integration. The parameter c is the electric permittivity of 
the medium. 


The relationship between the electric field E and the electric 
potential difference V ab between points b and a is given by the line 


integral 
< 36 > V ab 


*b 

E • dr , 

Ja 
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which is in accord with Eq. (30). The magnitude of the potential 
difference between two electrodes is usually called the applied 
voltage. 

There can be space-charge contributions to the field as well as 
surface-charge contributions, but for diffusing defect densities below 
10 vcm, for example, the space -charge contributions may be neglected 
for oxide thicknesses below 1000 A or so-? Because the space charge is 
a complexity in itself, in the present development we choose to neglect 
it. The present development is therefore restricted to the zero space - 
charge limit, so that the electric field is due to surface charge only. 
In the case of planar oxides, the zero space-charge limit leads to 
electric fields which are independent of position x within the oxide, 
and for this reason it has been designated^ the "homogeneous - f ie Id 
approximation" . At the moment, we restrict our consideration to the 
zero space -charge limit in order to focus more intensely on the purely 
geometric effects which is the subject of this presentation. 

B . Spherical Geometry 

Consider specifically the spherical geometry case with concentric 
spherical electrodes of radius a and b, with a < b. The electric field 
does not require the existence of electrodes, but it is helpful to 
imagine that they are present. Radial symmetry and zero space charge 
in the uniform medium leads to a ready evaluation of Gauss' law for 
this case, 

(37) £ 4 7r r 2 = , (a < r < b) 

where Q a is the charge on the inner electrode a. The radial component 
of the field is the only nonzero component, as can be argued 
convincingly from the symmetry of the physical problem. 


The corresponding electric potential difference evaluated at some 
arbitrary position r between the electrodes is given by 

r r Q' 

E dr «= r 1 
, r 4 tt£ 

Ja 

The total built-in electric potential across the oxide is given by 
this expression evaluated at r = b, 


(38) V(r) - V(a) 


f 1 - A 1 

1 r *J 


1 

b 


1 

a 


(39) V built _ in - V(b) - V(a) - ^ 

With no space charge, the charge Q at b is the negative of Q , 

d a 

(4°) Q b - - Q a . 

The specific evaluation of V built _ in will be carried out subsequently. 


C . Cylindrical Geometry 

Consider concentric cylindrical electrodes of radius a and b, 
with a < b. Gauss' law leads to 


(41) c E 2 7T r h = Q (a<r<b) 

r a 
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where the total charge Q on electrode a is determined by the charge X 
per unit length of that electrode, 

(42) Q a = X h . 

Conservation of charge and the absence of space charge allows writing 


(43) Q b - ■ Q a • 

As in the spherical geometry case, the radial component is the only 
nonzero component of the field. 


The electric potential 
is given by 


r 


V(r) for this cylindrically- symmetric case 


(44) V(r) - V(a) - - 


Ja 


E dr 
r 


-(X /2jre)ln(r/b) . 
3 - 


The total built-in electric potential across the oxide is given by 
this expression evaluated at r = b, 

(45) V built _ in - V(b) - V(a) - - (X/2 tt £ ) ln(b/a) . 

D. Planar Symmetry 

For completeness, let us give the analogous results for planar 
geometry. For the metal -oxide interface at a and the oxide -oxygen 
interface at b, the oxide thickness is given by L=b-a. Planar 
symmetry and zero space charge leads to a ready evaluation of Gauss' 
law, 


(46) e E x A = Q a , 

where r - Q /A is the charge per unit area. The corresponding 
electric potential is given by 

(47) V(x) - V(a) - (r/ £ )(x - a) . 

Evaluation of this expression at x - b gives the potential difference 
across the oxide. 


E. More Formal Treatment 

A more formal treatment can be based on Poisson's equation 

(48) -V 2 V = p/z , 

where p is the local volume charge density and £ is the electric 
permittivity of the oxide. For zero space charge, p=0. Employing 
E = -W, this leads to 

(49) -V 2 V --V-W-V-E-0. 

In spherical polar coordinates this becomes 
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(50) 


sin 9 


[ sin 6 fc (r2 V + r h { (sln 6 V + 4^] “ 0 


while in cylindrical polar coordinates this relation is 


(51) r f^ (r + - ^ + 


1 ^ 

r 39 


3 E z 
<9z 


Due to the spherical and cylindrical symmetry assumed in the 
respective cases, only the radial component of the field is nonzero. 
Furthermore, E - E (r) only, since, with the specified symmetries, 
there is no cause for variation of the field in the other directions. 
Thus we obtain in the spherically- symmetric case 

(52) |^(r 2 E r ) - 0 , 


while for the cylindrically- symmetric case we obtain 

(53) |^(r E r ) - 0 . 

These expressions are readily integrated. Thus for spherical 
symmetry, 

(54) E a , 

r r z 

while for cylindrical symmetry, 


(55) E oc — 
r r 


in accordance with the above results obtained using Gauss' law. The 
planar case is just as simple, since for zero space charge, 


(56) V 


E 


<3E V 3E V 3E- 
teT + 3y + dz~ 


0 . 


For symmetry over planes perpendicular to x, we expect that Ey-0 and 
E z =0, and that E x will be independent of y and z. Thus 



and 


(58) E = Constant . 
x 

The electrical potentials follow from a line integral of the 
electric field in the usual way, as already shown above. 


VI. Coupled Currents for Spherical and Cylindrical Symmetries 
A. Coupled- Currents Condition 

Let us consider the interface between the material being oxidized 
(usually a metal) and the oxide layer to be located by the radius a. 
This is the metal-oxide interface. Then the opposite interface of the 
oxide, which is the oxide-oxygen interface in contact with the 
attacking oxygen phase, is considered to be located by the radius b. 

As the oxide film grows, both a and b change. In spherical, 
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cylindrical, and planar geometries, the oxide film thickness L is given 
by the difference between a and b, so a = a(t) and b = b(t) . 

The electrochemical potential difference between the metal-oxide 
and the oxide -oxygen interfaces provides the driving force for the 
currents, as outlined above. At any position in the oxide layer, the 
existing electric field will aid transport of the ionic or the 
electronic species, and oppose transport of the other. The currents 
can themselves lead to a net charge transport, which changes the field 
since Q a will be modified. Careful analysis^ has shown that the net 
charge transported through the oxide layer of thickness L over any time 
increment is generally much smaller than the charge transported by 
either the ions or the electrons individually. Were this not the case, 
the electric fields produced within the system would reach catastrophic 
values very quickly, with dielectric breakdown occurring. The 
fundamental point is that a monolayer of charged particles can produce 
an incredibly large electric field, so that the net charge setting up 
the surface -charge field requires only a tiny imbalance in positive and 
negative charge. The difference between the positive and negative 
charge transport during the building of a layer of oxide can represent 
at most only a tiny fraction of the total particles in that layer. 

The coupled- currents condition 

(59) q x J l + q 2 - 0 

represents a state of zero net charge transport, since the charge 
current q^J^ of the ionic species (species 1) is nullified by the 
charge current q2 J 2 °f t ^ xe electronic species (species 2). This 
presumes the specific situation where there is one dominant diffusing 
ionic species and one dominant diffusing electronic species, which is 
the simplest case for growth by charged particles. 

The coupled currents condition has been used widely in the 
evaluation of planar oxide growth-? but to date it has not been applied 
to the growth of non-planar oxides. For the spherical and the 
cylindrical symmetries, the currents J (r) given by Eqs . (4) and (8) 
respectively, can be substituted into fq. (59) above to obtain the 
following form for the coupled-currents condition, 

(60) q 1 I 1 + q 2 I 2 = 0 . 

This is the most useful form for present purposes. 

B. Differential Equation for the Current 

The current density J s of species s is given in general, by Eq . 
(29). The geometry in question determines the form of the gradient 
operator V and the functional form of the electric field E. In both the 
spherical and the cylindrical geometries, C s varies only radially, so 
that C s = C s (r). The gradient operation in both geometries [see Eqs. 
(26) and (27)] then reduces to the especially simple form 

(61) VC g = r(dC s /dr) , 
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where r is the unit vector in the radial direction. In both 
geometries the electric field is radial, so that 

(62) E — r E r . 

Substituting these two relations into Eq. (29) yields the current 

(63) J = r J 

s s 

which is totally radial, with J g given by 

(64) J =-k n TB^s + qBEC. 

v 's B s d r s s r s 


VII. Solution of Problem for Spherical Symmetry 

A. Concentration Profile and Steady-State Current 

The electric field for spherical symmetry given by Eq. (37) can be 
substituted into Eq. (64) for the current, and in addition Eq. (4) can 
be used to replace the position- dependent current density J g by the 
position- independent quantity I . Thus we obtain 


(65) 


4 7rr ^ 


k T B 
B s dr 


+ q B C 


s s s 


Rearranging this equation to separate variables yields the integral 
form 


( 66 ) 


C s (r) 


'C s (a) 


(4irek B T^)dC s 
(*?s *%Qa)^s * c ^s 


dr 


These integrals are readily evaluated to give the concentration profile, 

(67) C g (r) - (cl s /q s % Q a ) + [ C s (a) - (cl s /q s % Q a ) ] P s (r) , 
where 

(68) P g (r) - exp [ a g [ (1/a) - (1/r) ] ] , 
with 


(69) a s = q s Q a /4^ £ k B T . 

The total particle current I s for species s then follows by evaluating 
the concentration profile expression at r«b, 

C s (b) - C s (a)P s (b) 


(70) I s - (, s Q a » A) 


1 - P,(b) 


Thus I is given in terms of Q and otherwise known parameters of the 
system for any species s. 

The subsequent algebraic details are simplified by introducing the 
parameter T defined by 
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(71) T - exp [ [ (1/ a) - (1/b) ] ] , 

where 


(72) a o = eQ a /4»rck B T , 

with e representing the electronic charge magnitude. Expressing the 
charge q s per particle of species s as 


(73) q g = Z s e , 

the quantity P s (b) appearing in Eq. (70) becomes 

(74) P g (b) - r Zs , 


since 
(75) a s 


Z a 
s o 


B. Built-in Potential 

To determine Q a , we now invoke the coupled currents condition. 
Substituting Eq. (70) for I s (with s-1 and s=2) into Eq. (60) and 
utilizing Eqs. (73) -(74) yields 


(76) 


Z 1 ®L t c i< a > rZl - C i< b >) 
r z l - l 


z2 


\ [C 2 (a)F 


Z 2 . 


C 2 (b) 


= 0 


- 1 


This represents an algebraic equation for T, and the solution yields 
the quantity Q a because that is the only variable quantity in the 
definition of r, as can be noted from Eqs. (71) and (72). This in turn 
is sufficient to evaluate V built _ in by means of Eq. (39). Thus 

(77) V built _ in - -(k B T/e)ln T . 

The details of solving the algebraic equation for T hinges upon 
the ratio Z 1 /Z 2 , but this is a technical problem only. The important 
point is that the solution will yield a value for T which depends upon 
the fixed parameters of the system but which under normal circumstances 
will be independent of the values of the radii a and b. Thus T will 
not depend upon the thickness of the oxide, and hence the built-in 
potential V built _ in will be fixed during growth of the oxide. This is 
a very important result for the kinetics of oxidation for the case of 
spherical geometry. 

At this point it is worthwhile to find the explicit forms for the 
built-in potential for some specific situations. Consider first the 
ionic species to be monovalent cation interstitials, in which case Z\ = 
1, with the attendant electronic species being electrons, so Z 2 - -1. 
The solution of Eq. (76) is readily obtained, 

(78) T- t » L C 1 (b) + S^C 2 (a)] / [J^CjU) + ^C 2 (b)] , (2^+1; Z 2 =-l) 

and the built-in potential follows immediately from Eq. (77). 
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Since this result is independent of the radii a and b, it follows 
that it must also give the corresponding value for the built-in 
potential in the planar limit where L~(b-a) is small relative to a. 
[Refer to ch. 7, §3 in Ref. 2 for the planar treatment]. 

For the converse case of monovalent cation vacancies and electron 
holes, « -1 and - +1. The same form of equation results, but 
with T replaced by l/T. Thus the solution for T for this case is the 
reciprocal of the right-hand side of Eq. (78) above. The built-in 
potential will then have the same magnitude but will be of opposite 
sign. 


Whenever the ionic species is divalent, the ratio Z^/Z^ will have 
the value -2. This case may involve cation interstitials (Z^*+2) with 
electrons, or else cation vacancies (Z^--2) with electron holes. In 
either case a quadratic equation for T emerges, but this is readily 
solved by the quadratic formula. The value of T will be independent of 
the thickness of the oxide, so the built-in potential given by Eq. (77) 
will be independent of oxide thickness. 

C . Ionic Current for Spherical Samples 

The total ionic current I]_ is readily evaluated at this point, 
since the evaluation of T is sufficient to yield P s (b) from Eq. (74) 
above, and also sufficient to yield a Q and hence Q a from Eqs . (71)- 
(72), 

(79) Q a = (47T£k B T/e) - £ ) l In T . 

Thus the charge on the metal changes with the radii a and b as the 
oxide grows. This technically is a violation of the coupled-currents 
condition, but it can be shown to provide only a minor perturbation on 
the final result. 

Utilizing the above result for Q a ,we obtain from Eq, (70) for 1^ 
the following expression for the total ionic current, 

(80) I x = k«Z x 3^ k B T ln(D [\ - g l 1 

Alternatively, the total current may be expressed in terms of the 
constant built-in potential, 

II ' 1 

(81) I x = -4^ ^ (- - g ) V built _ in x 

Cj(b) - C 1 (a)exp[-Z 1 eV built . ln /k B T] 

1 - exp[-Z ie V built . in /k B T] 

For the positive monovalent ionic species case, the evaluation of 
V b u i 1 1 -in g iven by Eqs. (77) -(78) can be inserted directly into this 
expression to obtain an alternate form, if desired. 

D. Oxide Growth and Particle Burnup Rate Laws for Spherical Samples 

Because 1^ gives the total current flowing through any symmetry 
sphere, as defined in Eq . (4), and moreover is independent of the 


C X (b) - C x (a)T 1 

1 . pZl 


XVI I- 15 



radius of the sphere chosen, it represents the value of the integral 
appearing in Eq. (3) for the volume rate of oxide formation. Thus we 
can write 


(82) 


dV, 


oxide 


dt 


R i h 


where is the total ionic current as evaluated above. Since 1^ is 
given in terms of the time -dependent radii a and b, it is necessary to 
relate these quantities to V oxide before the growth law can be 
evaluated. 


(83) 


Let us next define the volume expansion parameter cr, 
Volume of oxide formed 


Volume of metal used in oxide formation 

Since the volume of metal which has been utilized at any time t is 
given by 

4 3 4 3 

i 

o 


(84) V, 


metal 3 ^ a 


— 7T a 


where a Q is the radius of the metal sphere at t=0, the corresponding 
volume of oxide which has been formed is given by 


(85) 


'oxide 


- ( 


7 r b~ 


4 

— 7T a 


3 i (4 ,3 

J ' U * 


r 7T a 
3 o 


) . 


where b Q locates the oxide-oxygen interface at t=0. If no initial 
oxide is present, then b Q = a Q , and the above expression simplifies 
accordingly. The volume expansion parameter thus can be written as 

(b 3 - bo) + (a^ - a 3 ) 3333 

(86) a f f - 1 + [0> - bV(a^ - a J ) ] , 

a g . o o 

which provides the following relation between b and a, 


(87) b = (b 3 + (a - l)(a 3 - a 3 ) ] 1/3 . 


Whenever the volume of oxide exceeds the volume of metal consumed in 
forming the oxide, a exceeds unity. For a negligible expansion in 
volume upon oxide formation, a would be unity and b would maintain its 
t=0 value b Q . Since in general 

(88) V D x i j e — <7 V m e t a \ > 


we obtain from Eq. (84), 


(89) 


dV 


oxide 


dV, 


metal 


, 2 da 

-bnoa. — 
dt 


dt dt 

Setting this equal to Rill in accordance with Eq. (82) above gives 


(90) 


-era 


2 da 


^ f I _ I } = q 

l a h ) * 


dt ^ a b 

where Q is the constant quantity 

C x (b) - C 1 (a)exp[-Z 1 eV buiLt . in /k B T] 
1 - exp[-Z 1 eV built . in /k B T] 


(91) 0 


^built-in 
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with ji- being the ionic mobility defined by Eq. (33). Next we can 
utilize Eq. (87) for b in terms of a, and integrate the above 
differential equation from t = 0 to arbitrary time t, 

(92) | a [a 2 - a 2 ] + a ^(a) =0 t , 


where 


(93) 



da 


b Q + (ff - 1) (a Q - 


3 a 3 )F> 


The integral S (a) is easily evaluated, 

(94) ^(a) = - n J-_ 1} - [ (b^ + (a - l)(a 3 - a 3 ) ) 2/3 - b 2 ] 

We note in passing that 5 (a) can likewise be written in terms of b, 

< 95) \ ~ - 2(» X - 1) ' b o ] ' 

Combining the integrated form for S (a) with Eq . (92) yields the rate 

law for the time-rate of change of the radius a(t) of the metal 
particle . 

Next let us deduce the corresponding rate law for the variation of 
the outer radius b with time. From Eq. (87) which relates b to a we 
obtain by differentiation 

(96) a 2 da = -(a - l)' 1 b 2 db . 


Substituting this relation into Eq. (90) yields 
r a 1 T b db . db 1 _ 

<v> (;-n Hrs - b d? J- n ■ 

Integrating this equation from t-0 to arbitrary time t yields the rate 
law for the outer radius b(t), 

(98) )[^(b> - f [b 2 - b 3 ]]-Ot , 


where, employing Eq. (87), 
(99) S 2 (b) = £ 


b 2 db 


(a 3 - (a - 1) 1 [b 3 - b 3 ] ) 


1/3 


This integral is readily evaluated to obtain 

(100) * 2 (b) = - \ (a - 1) [ ( a 3 - (a - l)' 1 (b 3 - b 3 ) ) 2/3 - a 2 ] 
Note that * (b) can likewise be written in terms of a, 

1 r 2 2] 

(101) S 2 = - ^ (a - D [a - a Q J 

Combining the integrated form for -^>(b) with Eq. (98) yields the rate 
law for the time -rate of change of the radius b(t) of the metal 
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particle. This result also can be obtained by substituting directly 
Eq. (87) into Eq. (92) and utilizing Eq. (95). 


Next let us ask whether or not it is possible to obtain any 
expressions for the time -dependence of the oxide thickness L and the 
volume of oxide. Because L-b-a, the above expressions relating b and a 
Co t provide a method to obtain L(t) numerically. More directly, 
combining Eqs. (84), (85), and ( 88 ), the obtained relation 


(102) V oxide 






serves to yield V oxide directly in terms of a and also provides a 
ready way to obtain L in terms of a. Substituting b=L+a into this 
relation and solving for L yields 

(103) L = ( (1 - a) (a 3 - a*) + b^ ) 1/3 - a . 

Numerically the problem can then be evaluated as follows . A sequence 
of values of a is substituted into Eqs. (92) and (94) to obtain the 
corresponding sequence of time values . Next the corresponding 
sequences of values for b, V oxide , and L are obtained from the sequence 
of values of a by using Eqs. (87), (102), and (103). The results can 
then be plotted as individual curves, namely, a(t) vs t, b(t) vs t, 
V oxide (t) vs t, and/or L(t) vs t, as desired. 


There will eventually be a point in time where a(t) decreases to 
zero. Since the metal particle will then have been converted totally 
to oxide, we designate this as the burn-up time t burnup . The a(t) ■+ 0 
limit of Eqs. (92) and (94) yields the time required for the complete 
oxidation of the spherical particles , 


(104) t burnu p 


(a/ 2 0 ) 


2 


a 

o 



+ 


(a 




This can be written in an alternative form involving the quantity b roax 
defined as the value b approaches as a -* 0 . Then the total quantity 
of oxide is given by 


(105) 

4 K 3 

3 


so 

(106) 

ax “ 






) 



3^/3 

VJ 


Thus Eq. (104) can be expressed in the form 

(107) t burnup - (a/20) [a \ !)(*>«« ' h l )] ■ 

Due to the complexity of the above rate laws, little can be 
deduced from a cursory inspection. Three methods suggest themselves 
for obtaining physical understanding of the result, viz. , numerical 
evaluations using the scheme outlined above, series expansions, and 
simplifying approximations. Although we have used all three, at the 
moment it is most informative to adopt the latter approach. One 
simplifying assumption which can be made is that the volume expansion 
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of the oxide is relatively small. In the limit that no expansion of 
the metal occurs upon formation of the oxide, the volume expansion 
parameter a will be unity. The differential equation for a(t) given 
by Eq. (90) then reduces to 


(108) 


-a 


da r 1 1 i 

dt 1 fl " h. J 


0 


since the location of the oxide surface b(t) will then be unchanged 
with the degree of oxidation, having the initial value b 0 . 

Integration of this expression gives 

(109) [a 2 - a 2 ] - (2/3b Q )[ao - a 3 ] - 2 Q t . 

Since a Q < a for spherical particles, and b Q > 0, the second bracketed 
term tends to cancel the first bracketed term on the left-hand side of 
this equation. Furthermore, the value of the second term will always 
be smaller than the first term, so that the expression is well-behaved 
as time increases. The observation to be made is that the time 
required to reach any given a(t) will be shorter due to the presence of 
the second bracketed term. As shown later, the first bracketed term 
expresses the entire planar- limit result for which the radii a and b 
are much greater that the oxide thickness L-b-a. 


The simplified equation above yields the following result for the 
burn-up time for spherical metal particles, 

(HO) t burnup - a 2 /6fl , 

where for simplicity we have now assumed there exists no initial oxide 

(viz . , at t=0 , b - a ) . 

o o 

It is interesting to derive the corresponding result for the oxide 
thickness L(t) as a function of the oxidation time t. Again assuming a 
= 1 , we have 

(111) b = b Q = a + L , 
so that 

(112) da - -dL . 

Equation (108) takes the form 

(113) [1 - (L/b )] L dL = O dt , 

which is readily integrated to obtain the oxide growth rate law for 
the oxidation of spherical particles , 



Often it can be assumed that L 0 - 0. The cubic term is L is always 
smaller than the quadratic term, since L < b c , so this equation 
represents a well-behaved rate law. Furthermore, the cubic term 
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subtracts from the quadratic term, so a shorter time t is required to 
reach a given oxide thickness L because of the presence of the cubic 
term. 


It is interesting to compare this result with the corresponding 
result for planar geometry. The planar limit is achieved by assuming 
that L « b Q , so that (L/b 0 ) can be neglected relative to unity. The 
differential equation for growth [Eq. (113)] then reduces to the 
approximate form 

(115) L dL - O dt , 

which is integrated readily to yield the planar growth law, 

(116) L 2 - L 2 - 2 0 t . 

Comparison of this planar result with Eq. (114) for spherical samples 
shows that the cubic term contains all of the effects of the spherical 
geometry . 

If a metal plate having an initial thickness 2a 0 is attacked by 
oxygen from both sides , then a comparison between the planar result and 
Eq. (114) above for spherical samples shows that the burn-up time 
for bare spherical particles of initial radius a Q is only 1/3 the burn- 
up time required for the plate. Moreover, if the plate is attacked on 
only one side, the time required for the oxide to reach a thickness on 
the plate which is equal to the spherical particle diameter would be a 
factor of 12 larger than the burn-up time of the spherical particle. 
This simplified example thus serves to illustrate that geometrical 
effects on the rate of metal oxidation can be quite pronounced. 

A numerical computation has been carried out to illustrate the 
functional form of the simplified rate law [Eq. (114)]. The numerical 
values used for the several parameters are listed in Table 1. The 
results of the computation for oxide film thickness L vs time t are 
shown in Fig. 1. For comparison, the planar result given by Eq. (116) 
is also plotted in Fig. 1. The faster burn-up of the spherical 
particles relative to the planar sample as the oxide thickness 
approaches the value of the initial radius a Q of the metal particle is 
to be noted in the figure, corresponding to the factor of 3 deduced 
analytically above. In addition, the rapid increase in the rate as the 
oxide thickness approaches the value of the initial radius a 0 of the 
metal particle, corresponding to the burn-up point where the radius of 
the spherical metal particle approaches zero, is of major significance. 
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TABLE 1 


Values used for the parameters in computing Fig. 1. 


- 8 2 

= - 1 X 10 cm /V-sec 
- 6 2 

= 1 x 10 cm /V-sec 
14 3 

C 1 (a) = 1 x 10 /cm 

(^(b) = 5 x 10 17 /cm 3 

C 2 (a) = 1 x 10 17 /cm 3 

C 2 (b) = 5 x 10 17 /cm 3 

-22 3 

R = -2.0 X 10 cm 
T = 900 K 



a - 1 X 10" 4 cm =10,000 A 
o 

b = 1 x lo ' 4 cm =10,000 A 

O 

a — 1 

k„ - 1.3806 x 10‘ 23 J/K 
B 

e = 1.6022 x 10' 19 C 


Computed Results: 
V b uilt-in - -0.1210 Volts 
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Fig- 1. Oxidation rate law for spherical particles (upper curve) 
with growth by anion interstitial or cation vacancy 
diffusion coupled with electron-hole diffusion. The 
corresponding parabolic law [L 2 a t] for planar samples 
(lower curve) is shown for comparison. 
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